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Abstract. For inhomogeneous classical Coulomb fluids in thermal equilibrium, like the jellium or the two- 
component Coulomb gas, there exists a variety of exact sum rules which relate the particle one-body and 
two-body densities. The necessary condition for these sum rules is that the Coulomb fluid possesses good 
screening properties, i.e. the particle correlation functions or the averaged charge inhomogeneity, say close 
to a wall, exhibit a short-range (usually exponential) decay. In this work, we study equilibrium statistical 
mechanics of an electric double layer with counter-ions only, i.e. a globally neutral system of equally charged 
point-like particles in the vicinity of a plain hard wall carrying a fixed uniform surface charge density of 
opposite sign. At large distances from the wall, the one-body and two-body counter-ion densities go to 
zero slowly according to the inverse-power law. In spite of the absence of screening, all known sum rules 
are shown to hold for two exactly solvable cases of the present system: in the weak-coupling Poisson- 
Boltzmann limit (in any spatial dimension larger than one) and at a special free-fermion coupling constant 
in two dimensions. This fact indicates an extended validity of the sum rules and provides a consistency 
check for reasonable theoretical approaches. 
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PACS. 82.70.-y Disperse systems; complex fluids - 61.20.Qg Structure of associated liquids: electrolytes, 
molten salts, etc. - 82.45.-h Electrochemistry and electrophoresis 



1 Introduction 

One of relevant problems in soft condensed matter is the 
study of thermodynamic properties of charged mesoscopic 
objects (colloids or macro- ions of several thousand elemen- 
tary charges), immersed in a polar solvent such as water 
containing mobile micro-ions of low valence. The charged 
object together with the surrounding micro-ions form a 
neutral entity, known as the electric double layer, which is 
of intense theoretical interest, see reviews [Tl[21[31i3j . Elec- 
tric double layers are important also for predicting an ef- 
fective interaction between macro-ions in a solvent and, in 
particular, an "anomalous" like-charge attraction [5,164. 

The general classical system of mobile micro-ions can 
be formulated in two ways. The case "counter-ions only" 
corresponds to one species of equally charged micro-ions 
which neutralize the macro-ion charge. This is the situa- 
tion when a mesoscopic object, dissolved in a polar sol- 
vent, acquires an electric charge through the dissociation 
of functional surface groups, releasing in this way counter- 
ions into the solvent 0. In the case "added electrolyte", 
there is an infinite reservoir of ± charge pairs. 

In this paper, we restrict ourselves to a single electric 
double layer with counter-ions only, see the geometry in 
Fig. [T] The system is {/-dimensional {v > 2), inhomoge- 
neous along the x-axis and translationally invariant along 
all remaining v — l spatial axes. The large macro-ion is ap- 
proximated by a hard wall filled by a material of dielectric 
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Fig. 1. The general electric double layer with counter-ions 
only. For simplicity, we consider the homogeneous dielectric 
case with e' = e = 1 and the valence of counter- ions q = 1. 



constant e' in the half-space a; < 0. Its surface at x = 
is charged by a fixed uniform surface-charge density ere, e 
being the elementary charge and say cr > 0. The wall is 
impenetrable to the g-valent counter-ions moving freely in 
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the coiTLplementary half-space a; > 0. To simplify the nota- 
tion, we set the valence q — 1 ; to restore the g-dependence 
we have to substitute a -^ a/q and e ^- ge in final for- 
mulas. The counter-ions are classical point-like particles 
immersed in a solution of dielectric constant e. For sim- 
plicity, we consider the homogeneous dielectric case with 
e' = e = 1 (vacuum in Gauss units), without electro- 
static image forces. The system is in thermal equilibrium 
at temperature T, or inverse temperature /3 = l/(fcBT). 
Although the dimension three (3D) is of primary phys- 
ical interest, many works deal with the two-dimensional 
(2D) version of Coulomb systems interacting logarithmi- 
cally which are sometimes exactly solvable, at a specific 
free-fermion coupling [SlI^ ITUllTTlfT^ (for review, see [13]) 
or, in the case of the Coulomb gas of point ± charges, in 
the whole stability interval of temperatures [H1IT5] . 

The model became, due to its relative simplicity, the 
cornerstone for developing theoretical methods to general 
Coulomb systems. The weak-coupling (high-temperature) 
limit is described by the Poisson-Boltzmann (PB) mean- 
field approach [16 and by its systematic improvements 
via the loop expansion pniTSllT^I^] . The opposite strong- 
coupling (low-temperature) limit was pioneered by Rouz- 
ina and Bloomfield [3P and developed subsequently in [HI 
[23] ; a basic ingredient was the formation of a 2D Wigner 
layer of counter-ions close to the charged wall. In the field- 
theoretical approach of Netz and collaborators [MlHSHSl 
[27] . the leading strong-coupling behavior stems from a 
single-particle theory and the next correction orders cor- 
respond to a virial fluid expansion in inverse powers of 
the coupling constant. A comparison with Monte Carlo 
simulations shows the adequacy of the method to capture 
the leading strong-coupling behavior, but its failure for 
the first correction. Recently [28,29,30 , a strong-coupling 
theory starting from the 2D Wigner crystal was proposed 
which reproduces the leading single-particle picture of the 
virial method and at the same time gives the first cor- 
rection which is in excellent agreement with Monte Carlo 
data. 

For dense homogeneous and inhomogeneous Coulomb 
fiuids in thermal equilibrium, like the jellium or the two- 
component Coulomb gas, there exists a variety of exact 
sum rules which relate the particle one-body and two- 
body densities, for an old review see J31J. The necessary 
conditions for these sum rules is that the Coulomb fluid 
exhibits good screening properties, i.e. the particle cor- 
relation functions or the averaged charge inhomogeneity, 
say close to a wall, exhibit a short-range (usually expo- 
nential) decay to zero. The bulk charge-charge correlation 
functions satisfy, in any dimension, the zeroth-moment 
and second- moment Stillinger-Lovett conditions [321155] . 
In 2D, also the compressibility [M1I55] and higher-moment 
[5S1I57] sum rules are available explicitly. For semi-infinite 
domains constrained by planar wall surfaces, the density of 
particles at the wall is related to the bulk pressure via the 
contact theorem [38,39,40,41,42 . The Carnie and Chan where r = |r| and tq is a length scale. Dimension one has 
generalization to inhomogeneous fiuids of the second-momentspecial features and is excluded from the analysis. The def- 
Stillinger-Lovett bulk condition gives the dipole sum rule inition of the Coulomb potential ([T]) implies the character- 
[131133]. A sum rule for the variation of the particle den- istic small-fc singularity v{k) = 1/fc^ in the Fourier space. 



sity with respect to the external electrical field was de- 
rived by Blum et al. g5]. The WLMB (Wertheim, Lovett, 
Mou, Buff) equations [^H7] . originally derived for neu- 
tral particles, were adapted to charged systems as well 
[35] . The charge-charge correlation functions decay slowly 
as an inverse-power law along the wall [1^150] and a sum 
rule for the amplitude function holds ^5r,52,53 . A rela- 
tion between this algebraic tail and the dipole moment 
was found in [51]. 

The aim of this work is to investigate whether the 
known sum rules apply to our electric double layer with 
counter-ions only. The problem is that the present model 
does not exhibit good screening properties: the density 
profile and the amplitude function of the long-range tail 
along the wall exhibit a long-range (inverse-power) tail at 
large distances from the wall. We check the sum rules in 
two exactly solvable cases: the PB theory in any dimension 
v > 2 and at the special free-fermion coupling /3e^ = 2 in 
2D. In both cases, all sum rules hold which indicates their 
extended applicability. The sum rules provide exact infor- 
mation about the model and represent strong constraints 
which have to be fulfilled within a reasonable theoretical 
description. 

The paper is organized as follows. The known sum 
rules are listed in sect. |2J The PB treatment in sect. |3] 
implies the inhomogeneous density profile; the correlation 
functions are evaluated with this profile by using a method 
which is alternative to the one-loop level of the field the- 
ory [^. The complete solution of the 2D model at the 
free-fermion coupling /3e^ = 2 is the subject of sect. |31 
Sect. [5] is the Conclusion. 



2 List of sum rules 

We consider the homogeneous dielectric case e' = e = 1 oi 
the electric double layer pictured in Fig. 1. Particles are 
constrained to the semi-infinite i/-dimensional Euclidean 
domain. Each point of this domain r = (a;, R) has the 
component x G [0,oo], along which the system is inho- 
mogeneous, and the remaining v ~ 1 perpendicular (un- 
bounded) components R S [—00,00]""^, along which the 
system is homogeneous and translationally invariant. 

In dimension i/, the Coulomb potential v at point r, 
induced by a unit charge at the origin 0, is defined as the 
solution of the Poisson equation 



Av{r) = -s^S{r), 



(1) 



= 2TT''/'^/r{iy/2) (in particular, S2 = 2tt and 
is the surface area of the z/-dimensional unit 
sphere. Explicitly, 



where s^, 
S3 = 47r) 



v{r) = 



-Hr/ro) iiv = 2, 
r'^~^ l(y — 2) otherwise. 



(2) 
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This maintains many generic properties, like screening and 
the corresponding sum rules, of "real" 3D Coulomb sys- 
tems with 1/r interactions. Dimensions v — 2,3 are of 
special physical interest. The 2D logarithmic potential cor- 
responds to the interaction of infinite uniformly charged 
lines which are perpendicular to the given plane. Since 
the point-like particles possess the same charge, no short- 
range regularization of the Coulomb potential is needed 
and the thermodynamics is well defined. 

The Hamiltonian of particles i = 1, . . . ,N with charge 
— e at positions {r-i = [xi^Hi)} and the fixed surface 
charge density ae sd, x ~ reads 



H 



i<j 



vi\r^-rJ\) + e a— 



Su v^ 



2 ^^* 



(3) 



In the canonical ensemble with the requirement of the 
overall neutrality, the thermal average at the inverse tem- 
perature /3 will be denoted by (•••). At one-particle level, 
we introduce the particle density 

n(r) ^ n{x) = (E.'5(r - r,)) . (4) 

At two-particle level, we have the two-body densities 

n2(r, r') = n2{x, x' ] |R - R'|) = 713(3;', x; |R - R'|) 

V^/(r-rO<5(r'-r,)). (5) 



The corresponding (truncated) Ursell functions 
712 {v,y') = n2{v,r') — n{v)n{r') 



(6) 



go to zero at large distances |r — r'| -^ 00. In the above 
formulas, we used invariance with respect to translations 
along the wall surface and rotations around the x axis. 

Now we list all known Coulomb sum rules, adapted 
to our system, forgetting for a while that the necessary 
conditions like screening do not apply. 



2.1 Sum rules involving the particle density only 

The overall electroneutrality of the system is equivalent 
to the condition 



dxn{x) = a. 



(7) 



Thus the particle density n(x) must go to zero faster than 
l/x as a; — ?► 00. 

The contact theorem for planar wall surfaces [5511551 
HDlHTUi^ relates the contact density of particles n(0) to 
the bulk pressure of the fluid P as follows 



/3P = n(0)--;3e2s,a2^ 



(8) 



Since the density of particles vanishes in the bulk x ^- cx), 
we have P = and the exact constraint 



n(0) = -/^e^s^CT^ 



(9) 



2.2 Electroneutrality and dipole sum rules 

Two sum rules for inhomogeneous Coulomb fluids have the 
origin directly in the zeroth-moment and second-moment 
Stillinger-Lovett conditions for the bulk correlation func- 
tions [5^155]. 

The electroneutrality condition with a particle fixed at 
some point takes the form 



n(x) 



f'"'!' 



dRn^2^\x,x'-R). 



(10) 



The dipole sum rule 



dx dx' dR{x'~x)n^p{x,x';R) 
Jo Jo J 



1 



(11) 

follows directly from the Carnie and Chan generalization 
to inhomogeneous fluids of the second-moment Stillinger- 
Lovett bulk condition [431I44J . Note that the integrals over 
X and x' in (|lip cannot be interchanged since the inte- 
grated function is not absolutely integrable; if the inter- 
changing of integrations would be possible, the result is 
zero since (x' — x) changes its sign under the exchange 
X O a;' transformation. 



2.3 Sum rule of Blum et al. 

Blum et al. |45) derived a sum rule which relates the vari- 
ation of the particle density n(r) with respect to the ex- 
ternal electrical field E — s^ae to the dipole moment with 
respect to a particle fixed at r. For our system, the sum 
rule reads explicitly as 



dn(x) 
da 



= -s,yPe 



I dx' f dR{x' ~x)n^^\x,x';R). 



(12) 

This relation is superior to the dipole sum rule dTTj) which 
results by integrating both sides of ((T2|) over x G [0, 00] 
and then considering the electroneutrality condition ([7|). 



2.4 The WLMB equations for charged systems 

The WLMB equations PSIHT] were originally derived for 
neutral particles. Their generalization to charged systems 
[5nH5] relates the gradient of the density to an integral 
of the corresponding two-body Ursell function over the 
boundary: 



dn{x) 
dx 



f dRn^P{0,x;R). 



(13) 



Note that integrating this formula over x £ [0, 00] and 
taking n{x — > 00) = 0, we recover the electroneutrality 
condition ([TU| for the special case x = 0. 
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2.5 Long-range decay along the wall 

Because of asymmetry of the screening cloud around a par- 
ticle sitting near the wall, the two-body Ursell functions 
decay slowly along the wall [i^lSU] . Using linear response 
in combination with a simple macroscopic argument based 
on the electrostatic method of images 51,52,53 , one ex- 
pects an asymptotic inverse-power law behavior 

n^^\x,x';R)^i^^, R ^ ^. (14) 

In the {v — l)-dimensional Fourier k-space with respect to 
R, defined by 

g{x,x';k) = / dRexp(— ik • R)g(a;,a;'; i?), 

/dk 
-———^exp{ik-'R)g{x,x';k), (15) 

this behavior is governed by the kink at k = [55 of the 
small wave number behavior 

h^P {x, x'; k) = n^P{x, x' ■ 0) - y /(x, x')k + • • • . (16) 

The function f{x,x'), which is symmetric in x and x' , 
obeys the sum rule [5n[52l[53] 



potential (f) at distance x from the wall satisfies the Poisson 
equation 

d^(i)(x) 



dx^ 



= -s^p{x). 



(20) 



In the mean-field approach, the particle energy at x is 
approximated by —e(f){x) in the corresponding Boltzmann 
factor: 

n{x) = nocxp [/3e(/)(x)] . (21) 



The boundary condition for the electric field reads 
dct>{x) 



dx 



— —SuCJe < 0. 



We introduce the Gouy-Chapmann length 

2 



Pe^aSn 



(22) 



(23) 



and the dimensionless electric potential ip{x) = /3e0(x), 
n(x) = no exp[V'(a;)], which satisfies the equations 



^^5=p[«.„, 



Their solution is 



dip(x) 



dx 



x=0 



ip{x) = c- 21n(a; + &), 



nge 



ab. 



l (24) 



(25) 



oo poo 



JQ 



dx dx' f{x,x') = - q 2 2 ' 



(17) The corresponding particle density 



It is interesting that for all exactly solvable cases the func- 
tion f(x,x') takes the product form 



n{x) = 



ah 



{x + hf 



(26) 



f{x,x') = -g{x)g{x'). 



(18) 



The decoupling of coordinates x and x' is intuitively due 
to the fact that the lateral distance between the points r 
and r' goes to infinity and so the x-coordinates of the two 
points become mutually uncorrelated. 

A relation between the algebraic tail of the correlation 
function along the wall and the dipole moment of that 
function was found in Ref. |54) : 

/•oo /• t^OC 

/ dx' / dR(x' — x)n2 {x^x';R) — ^ I dx f{x,x) 

for any x > 0. The integration of this relation over x £ 
[0,oo] leads to the equality which is consistent with the 
previous sum rules (|TT|l and (flT)) . 



evidently satisfies both the electroneutrality condition ([7]) 
and the contact theorem ^. Note that the spatial di- 
mension V enters only through the definition of the Gouy- 
Chapmann length (1231) ■ 

The correlation functions can be evaluated with the 
density profile (pS)) by using the one-loop level of the field- 
theoretical method of Netz and Orland [5D] . The correla- 
tion functions are introduced there as auxiliary quantities 
to deduce the first correction to the density profile. Their 
exact meaning is not specified. This is why we propose an 
alternative derivation. It follows from the general theory 
of Coulomb fluids, namely the Ornstein-Zernike equation 

h{T, r') = c(r, r') + ( dr"c(r, r")n(r")/i(r", r'), (27) 

which relates the truncated pair correlation 



3 Poisson- Boltzmann theory 

We adapt to arbitrary v >2 dimensions the derivation of 
the 3D PB theory by Andelman [T2] which was formulated 
originally longtime ago by Gouy |56j and Chapman [57j . 
The density of counter-ions n(x) corresponds to the 
charge density p{x) — — en(x). The average electrostatic 



/i(r,r') = 



n(r)n(r') 



(28) 



to the so-called direct correlation function c. For particles 
interacting via the Coulomb potential w, the leading term 
of c reads l50ll36l[37l 



c(r,r') ~ -Pe^v{y,v'). 



(29) 
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Inserting this expression into (I27p . applying to both sides 
the Laplacian with respect to r and using the Poisson 
equation ([T]), we obtain 

[Z\r - (3e'^s^n{x)] h{r, r') = (3e'^s^S{r - r'), x, x' > 

(30) 
and 

Arh{r, r') = 0, x <0,x' >0. (31) 

In the Fourier space with respect to R (|15p . we get 



After some algebra, the small-fc expansion of the truncated 
pair correlation is found to be 

h{x, x'; k) _ lx\+ ibx\ + ib'^x< + 36^ 
~ ~3" 



Pe^s, 



{x + h){x' + h) 

1,4 



{x + h){x' + b) 



(41) 



dx"^ 



{x + hf 



h{x, x' \ k) = Pe'^Si^6{x — x'), 



x,x' >0 (32) 



and 



dx^ 



h{x,x';k) =0, X <0,x' > 0. 



(33) 



where x< — min(a::,x'). 

The partial Fourier transform of the two-body den- 

sity ^2 {x,x']k) = n{x)n{x')h{x,x';k) has the small-fc 
expansion 

„(T), , _ 4 x^ -H Sbxl- + 3b'^x< + 3b^ 

"2 ^''''''^^""3/3e2s, (a; + 6)3(^/ + 5)3 

4 6^ 



These equations are supplemented by the conditions that 
h{x,x';k) and dh{x , x' ; k) / dx are continuous at a; = 
and by regularity h{x, x'; fc) — > as x — ;■ ±cxd. The formal 
solution of Eq. ^ reads [35] 



^;9e2s, (x-f fo)3(x' + 5)3^^'"' ^^^^ 
The first term of this expansion determines the integral 



f dIln^^\x,x';R) 



^^^^ = -![,.+ (x)^-(x')^(x-x') 

+f+{x')'f^{x)0{x' - x)] + Aip+{x)ip'^{x), 

(34) 

where A is a free parameter, 9 the Heaviside step function 



4 xl-+ 36x^ + 3b'^x< + 36^ 
3/3e2s^ (x-H6)3(x'-h6)3 ' 



(43) 



The presence of the second term in (|^^ signalizes the 
asymptotic behavior of type (J14p . With regard to formula 
p^ . the function /(x,x') has the long-range form 



/(x,x') 



(44) 



0(x) 



1 if X > 0, 
if X < 0, 



(35) 



the functions ip are given by 
d^ 2 



_dx2 (x-|-6)2 
and the Wronskian 



(p=^(x) =0, ^=^(±00) =0 (36) 



W = ^-{x)-^^+{x)~^+{x)^p-{x) (37) 

dx dx 

does not depend on x, as one can prove directly by differ- 
entiating W with respect to x and then using differential 
equations ([55)1 for ip"^. The solutions of Eqs. ([55]) read 



^±(x) 



1± 



1 



fc(x + &) 



exp [=Ffc(x + b)] (38) 



and the Wronskian W = — 2fc. The formal solution of Eq, 
dSSl) is 

h{x,x'; k) 



(3e^sl{x + b)^{x' + b)3' 

Notice that the product form ([T5| takes place with ^(x) = 
b'^y/8/{l3e'^sl)/{x -1-6)3. Using the last two relations it is 
easy to show that the sum rules (dU]), (HI]), (US]), (US]), (HZ]) 
and ([TO]) hold, in spite of the long-range nature of the 
density profile n{x) and the function /(x, x'). 

4 2D model at the free-fermion coupling 

The 2D version of the present model is mappable onto 
free fermions at the special coupling F = (ie^ = 2 [SJ 
150] . Here, we consider the Grassmann formulation of the 
Coulomb system on the surface of a semi-infinite cylinder 
of circumference W ^ i.e. the strip with x £ [0, cx)] and 
periodic boundary conditions along y £ [0,W] ^U\. The 
charge line density ae at x = is neutralized hy N = 
aW particles of charge — e. At the coupling F = 2, the 
renormalized one-body potential acting on particles reads 



Bip [x )exp(fcx). 



(39) 



Wr: 



.(x) 



An 



■ exp 



2-K ^ 



(45) 



where B is a free parameter. The continuity conditions at 
X = determine the parameters A and B as follows 



A = - 



1 



,2feb 



1 



In terms of two sets of anticommuting Grassmann vari- 
ables ^j and '\\)j, defined on a discrete chain of N sites 
j = 0,l,...,iV — 1, the partition function of the system is 
expressible as 



2fc l + 2(fc6)-f 2(fc6)2' 



B 



M 



{kbf 



k 1 + 2{kb) + 2{kby 



(40) 



■■^N 



Vil)V£, e 



S(?,V) 



^(C,^) 



N-l 



^^ ^jWjTpj, (46) 
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where the diagonal interaction matrix has the elements 



/"OO 

Wj = W / dx Wrcn{x) exp 
Jo 



second product Sjk'Sj'k produces the Ursell function of the 
form 



--jx]=^. (47) 



AT) 



{x, x']R) = —n 



x + x' + \R\ fx + x'-iR 



Due to the diagonalized form of the anticommuting action 
S'(^, ^), the partition function is available explicitly as the 
product of diagonal matrix elements 



Zjsi = WqWi ■ ■ ■ WN- 



(48) 



Within the Grassmann formalism, the particle density 
is given by 



N-l 



An 



W^ 



n{x) = w,en{x) Y^ {ijipj) exp ( --j^ja 
where the two-correlator 



(49) 



Zn J owj Wj 

(50) 
In the thermodynamic limit A^ — > oo, which is equivalent 
to VK ^- oo keeping the ratio N/W = a fixed, we pass 
from the semi-infinite cylinder to a 2D half-space. The 
variable t — {j + ^)/{N — i) becomes continuous and Eqs. 
|), (|50l) imply the density profile 

n{x) — Ana / dt t cxp(— 47rcrte) 
Jo 

r 

= An I ds s exp(— 47rsa:) 
Jo 



Anx"^ 



[l - (1 + Anax) e-^'T.Txj 



(51) 



Like in the PB approach, the density has a long-range 
tail: n{x) ~ l/(47ra;^) as a; — )■ cxd. It is easy to verify that 
that the electroneutrality condition ([7]) and the contact 
theorem ©, taken with /3e^ — 2 and S2 = 2n, hold. 
The two-body densities are given by 

N-l 
n2{r,r') = Wrcn{x)Wrcn{x) ^ (O'/'feG' V'fc') 



U+i' = k + k') 



2n, 



X exp [ - — [jz -f fcz -K i'z' + k'z'] ) , (52) 



where z — x + iy and z ~ x — iy are the complex coordi- 
nates of the point r — (x, y). Due to the diagonal form of 
the action S{^,ijj), the four-correlator {£,j'4'k£.j'''l'k') can be 
calculated using the Wick theorem, with the result 



valid universaly for the finite- M^ cylinder as well as for 
W ^>- oo 2D half-space. The integral representation of 
this formula for 2D half-space 



,(T)/ 



{x,x'-R)^~{Anf dtte~2"*("+"'M ds 



^^-2tts(x+x') ^-2TTm(t-s) 



(55) 



is useful to derive the important integral 

/OO n(7 

dRn^P{x,x'-R) = -{Anf / duV4-*(-+"'). 
-cxj Jo 

(56) 
Using this formula, the sum rules ([TOl), (HI]), dH]) and ([T3| 
follow immediately. The explicit representation of the re- 
lation ([S3]) for 2D half-space 



7r2[(x + x')2+i?2]2 

X |l - 2 [(1 + 2n(j{x + x')) cos.{2naR) 

+2naRs\x\{2naR)] exp[-27rcr(a; + x')] 
+ [(1 + 2ncf{x + x')f + {2n(TRf] 

X exp[-47rcr(x + x')] \ (57) 

is convenient to specify the asymptotic R ^- oo tail. We 
obtain the expected v — 2 behavior ((T4]) with the function 



f{x, x') = -4cr^ exp[-47ro-(a; + x')] 



(58) 



which is again of the product form (TTSl) with g{x) — 
2(jeyi\){—Anax). In contrast to the PB result (^1)) . this 
function is short-ranged with the exponential decay gov- 
erned by the surface charge density. The sum rules ([T7| 
and Hn]) are fulfilled. 



(OV'fcO'V'fc') 



1 



WjWj. 



{SjkSj'k' — Sjk'Sj'k) 



(53) 



The first product of Kronecker functions SjkSj'k' leads in 
([5^ to the term n{x)n{x') which has to be subtracted 
from 712 (r, r') to obtain the truncated Ursell function. The 



5 Conclusion 

The model of our present interest was the electric dou- 
ble layer with counter-ions only. Such system is "sparse" 
in the sense that the particle density and two-body den- 
sities vanish at asymptotically large distances from the 
wall. Moreover, the large-distance decay is usually not fast 
enough as required by the linear response argument in the 
derivation of standard sum rules for classical inhomoge- 
neous Coulomb fiuids. 

We studied two exactly solvable cases of the model. 
In the Poisson-Boltzmann limit, the long-range behavior 
takes place regardless of the dimension i/ > 2, for both the 
density profile (|26p and the asymptotic amplitude function 
(UU. As concerns the 2D version of the model at the free- 
fermion coupling F = 2, the density profile ([ST]) has the 
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long-range tail while the asymptotic function (|58p is short- 
ranged. 

The fact that the known sum rules are reproduced in 
two exactly solvable cases of the electric double layer in- 
dicates their extended applicability. The dipole sum rule 
PT|) deserves a special attention: the tendency of the sys- 
tem to its bulk regime with the second-moment Stillinger- 
Lovett condition was crucial in its derivation. As was men- 
tioned before, the bulk regime of our model is "emptiness" . 

The asymptotic amplitude functions (l44l) and (j58l) de- 
couple themselves in x and x' coordinates, in analogy with 
the jclhum and Coulomb-gas models [IHIEI]- 

We omitted from discussion "dense" systems like the 
one-component plasma with a neutralizing background 
(jellium) or the two-component Coulomb gas of ± charges. 
These systems exhibit exponential screening and so the 
sum rules automatically hold; their verification in special 
cases was given e.g. in [T^HSIIST] . 

The extended validity of general sum rules for inhomo- 
geneous Coulomb fluids can serve as a useful check for the 
adequacy of weak-coupling physical theories. The strong- 
coupling theories like ^511^150] are based on the Wigner 
ground-state lattice structure formed by counter-ions. The 
Coulomb system is not in its fluid phase in that regime 
and the sum rules do not longer apply. 

The support received from Grant VEGA No. 2/0049/12 is ac- 
knowledged. 
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